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Q. 

V«J ' Vibrations of the JefFcott rotor are modelled by a three degree of freedom system 

f^ . including coupling between lateral and torsional modes. The crack in a rotating 

shaft of the rotor is introduced via time dependent stiffness with off diagonal cou- 
plings. Applying the external torque to the system allows to observe the effect of 
crack "breathing" and gain insight into the system. It is manifested in the com- 
►^ , plex dynamic behaviour of the rotor in the region of internal resonance, showing a 

(^n I quasi-periodic motion or even non-periodic behaviour. In the present paper report, 

£2 ' 'we show the system response to the external torque excitation using nonlinear anal- 

og . ysis tools such as bifurcation diagram, phase portraits, Poincare maps and wavelet 

,jj; I power spectrum. In the region of resonance we study intermittent motions based on 

(^ • laminar phases interrupted by a series nonlinear beats. 
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1 Introduction 



In many technical devices and machines possessing rotors, i.e. turbines, pumps, 
compressors, etc., their actual dynamic condition determines their proper and 
safe operation. Often, these machines have to work over the extended pe- 
riods of time in various temperature regimes under large variable loading. 
As a consequence, the components of the machine are exposed to potential 
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damage, such as for example transverse crack [1,2]. For reliable and safe op- 
eration of such machines their rotors should be systematically monitored for 
presence of cracks and tested using non-destructive tests or vibration based 
techniques [21]. Structural health assessment of the rotating components may 
be performed by examining the dynamical response of the system subjected 
to external excitation applied to the rotating shaft. Such a strategy requires 
solution and study of the dynamical behaviour of the rotor [3,4]. The main dif- 
ference between the cracked and uncracked rotors can appear in the region of 
internal or combination resonances [5,6]. One should notice that the dynamics 
of the rotor is related to the coupling of torsional and lateral vibrations and 
the whole system is nonlinear. The crucial point is related to the model of the 
"breathing" crack, indicating that cracks involve additional parametric exci- 
tations with the frequency equal to an angular velocity of the rotor rotation 
[8,7]. The crack introduces extra nonlinearity terms and produces coupling 
terms of lateral and torsional modes. 




Fig. 1. Schematic plot of the cracked Jeffcott rotor [7]. 



Tondl [9], Cohen and Porat [10] and Bernasconi [11] showed that typical lat- 
eral excitations, such as unbalance, may result in both lateral and torsional 
responses of uncracked rotor. Continuing these works Muszynska et al. [12] 
and Bently et al. [13] discussed rotor coupled lateral and torsional vibrations 
due to unbalance. Their experimental results exhibited the existence of sig- 
nificant torsional vibrations, due to coupling with the lateral modes. Further 
reported works have used this coupling under external excitation as a means 
to identify the presence of cracks [14,15]. It has become clear that the spec- 
tral components of the shaft response under radial excitation can be used to 
identify the presence of transverse shaft cracks. In particular, Iwatsubo et al. 
[16] considered the vibrations of a slowly rotating shaft subject to either peri- 



odic or impulsive excitation. They identified specific harmonics in the response 
spectrum, which are combinations between the rotation speed and excitation 
frequency and they can be used to detect the presence of the crack. Also, Iwat- 
subo et al. [16] noted that the sensitivity of the response to the magnitude of 
the damage depends on the value of the excitation frequency chosen for the 
detection. Finally, Ishida and Inoue [17] considered the response of a horizon- 
tal rotor to harmonic external torques. Furthermore, Sawicki et al. [18,19,7] 
examined transient response of the cracked rotor, including the stalling effect 
under the constant driving torque. 

In this paper we continue these investigations. By applying the external har- 
monic torque excitation we show the characteristic features of internal reso- 
nances leading to quasi-periodic and non-periodic types of motion. 



2 The model and equations of motion 



The model of the Jeffcott rotor, which is by a definition limited to a rotating 
mass disc and massless rotating elastic shaft, is presented in Fig. 1. 

The angular position of the unbalance vector can be written as 



d{t)=ujt + ij{t)+(3, (1) 

where a; is a constant spin speed of the shaft, ■?/'(t) is a torsional angle, while 
j3 is the fixed angle between the unbalance vector and the centerline of the 
transverse shaft surface crack. It should be noted that 



^ = $ = cj + V', e = i), (2) 

where $ is the spin angle of the rotor (Fig. 1). 

The kinetic T and potential U energies for the rotor system subjected to lateral 
and torsional vibrations can be expressed as follows: 
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T = - Jp^2 ^ _^ (^^2 ^ ^2^ ^ -Mt^e"^ + MeO {-z sind + ij cos 9) (3) 

U = ]^{zy}}^A'\ + h,t^\ (4) 

where the stiffness matrix reads 



K, 



h h 

'^yz i^yy 



(5) 



The coupled equations for transverse and torsional motion for the rotor system 
take the following form: 



Mz + Ciz + k^^z + kzyy = F^ + Me U)"^ cos + sin 
My + CtZ + ky^z + kyyy = Fy + Me (i)"^ sin + cos 



(6) 

(7) 

Jp^ + Ct-ip + eCi{z sm0 -y cos 0) - e{F^ sin - Fy cos 0) (8) 

dU 

+e [{k^^z + k^yy) sin - {ky^z + kyyy) cos 6*] + — = Te, 

where M and Jp is the mass and mass moment of inertia of the disk, while e is 
eccentricity of the disk. F^ and Fy are the external forces (including gravity) 
in z and y directions, respectively, Ci and Ct are lateral and torsional damping 
coefficients, Te is an external torque. In our case 



-Mg, Fy = 0, Te = Tosin(cjet). 



(9) 



where g is the gravitational acceleration, Tq and uo^ are amplitude and fre- 
quency of the external torque, respectively. 

The stiffness matrix for a Jeffcott rotor with a cracked shaft in rotating coor- 
dinates can be written as: 



K 



R 



k^ 

U K/ff 




(10) 



where the first matrix, on the left hand side, refers to the stiffness of the 
uncracked shaft, and the second defines the variations in shaft stiffness k^ and 
kr, in ^ and 77 directions, respectively. The function /($) is a crack steering 
function which depends on the angular position of the crack $. The hinge 
model of the crack might be an appropriate representation for very small 
cracks, Mayes and Davies [20] proposed a model with a smooth transition 
between the opening and closing of the crack that is more adequate for larger 
cracks. In this case the crack steering function, or the Mayes function, takes 
the following form: 



fm 



1 + cos(<l>) 



'111 



The stiffness matrix for a Jeffcott rotor witli a cracked sliaft in inertial coor- 
dinates, K/ is derived as 



K, = TKrT-\ T 



cos $ — sin $ 
sin $ cos $ 



(12) 



wliere 



Aki 



Ak^ + A^ 



Aki 



Ak^ - A^ 



(13) 



Consequently after [7]: 



dU kdf{^) 



{Aki + Ak2Cos2^)z^ 



dtp 4 d^p 

+2zyAk2 sin 2$ + {Aki - Ak2 cos 2$)i/^ 

, f{^)kAk2 



z sin2<l> - 2zt/cos2<l> -y sin2<l> + hip 



(14) 



3 Results of Simulations 



Assuming the example of simply supported steel shaft, with a disc at the 
center and a transverse crack located near the disk, we applied Eqs. 6-8 and 
performed the simulations. The parameters of the examined system are shown 
in the Table 1. Simulations were done using the Euler integration procedure 
with a integration step 6t = 27r/(36000c<j) and a sampling time At = lOOO^t. 

The results for the range of rotational velocities uj G [100, 250] are presented 
in the bifurcation diagram (Fig. 2). One can see that the resonances appear 
as they are expected just below the lateral natural frequency and slightly 
above the half of lateral natural frequency, which one would expect for a 
non-linear system of this kind. Spanning the rotor spin speed in the range of 
145-149 rad/sec and 216-234 rad/sec, reveals the presence of a wide variety 
of the excited response characteristics. To explain the nature of vibrations in 
the regions of resonance we focused on the specific frequencies denoted by 
'l'-'9' for the corresponding values uj = 120, 132, 133, 140, 180, 193, 200, 
201, and 220 rad/s. In Fig. 3a we show the time histories of displacement 
in the vertical direction for nonresonat cases u = 120, 140, 220 rad/s. Black 
points are stroboscopic points plotted with the rotation frequency ui. Note 
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Fig. 2. Bifurcation diagram of the cracked rotor. The numbers 'l'-'9' correspond 
to w = 120, 132, 133, 140, 180, 193, 200, 201, and 220 rad/s. 
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Fig. 3. Time histories of displacement in the vertical direction (a) for to = 120, 140, 
220 rad/s (small amplitude periodic vibrations) and (b) for uj = 132, 200, 201 rad/s 
(intermittent vibrations). Black points shows stroboscopic points. Numbers '1', '4', 
'9', '2', '6', and '7' indicate the corresponding values in the bifurcation diagram 
(Fig. 2). 

that all three considered here cases indicate periodic vibrations with a small 
amplitude. On the other hand in Fig. 3b we show the resonant time series 
with rotation velocities u = 132, 200, 201 rad/s. In contrary to the previous 
cases the displacement z can be large. Interestingly it increases after small 
laminar behaviour with a relatively small amplitude. Note that in both panels 
the numbers '1', '4', '9', '2', '6', and '7' indicate the corresponding values u in 
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Fig. 4. Phase diagram projections into z-y plane for few values of speed velocities 
{uj = 120, 132, 133, 140, 180, 193, 200, and 201). 

the bifurcation diagram (Fig. 2). 

To shed some more hght into the lateral modes we plotted, in Figs. 4a-h 
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Fig. 5. Time histories of torsional angular displacements (j) for intermittent oscil- 
lations {io = 132, 193, 200, 201). Black points shows stroboscopic points. Numbers 
'6', '7', '8', '2' indicate the corresponding values in the bifurcation diagram (Fig. 1). 



the phase diagram projections into z-y plane are plotted for few values of 
speed velocities( lu = 120, 132, 133, 140, 180, 193, 200, and 201). Black points 
indicate the corresponding Poincare maps. In cases of nonresonant vibrations 
(Figs. 4a,d,e) of small amplitudes we observe with closed lines and a number 
of singular black points indicating periodic motion. Interestingly in Fig. 4c 
the Poincare map show a line instead of points indicating the quasiperiodic 
character of motion. The rest of plots (Fig. 4b,f,g) have non-periodic character, 
while Fig. 4h has features of moth quasi-periodic and non-periodic. Note that 
Fig. 4b,f,g, and h look like chaotic attractors, however to tell more about 
the non-periodic beats nature one should perform standard calculations of 
Lyapunov exponents or correlation dimension analysis. 



Finally in Fig. 5 we plot the time histories of torsional angular displacements ( 
u = 132, 193, 200, 201). Black points shows stroboscopic points. Numbers '6', 
'7', '8', '2' indicate the corresponding values in the bifurcation diagram (Fig. 
1). These plots shows the relaxation nature of the torsional oscillations. This 
is consistent with results for intermittent lateral oscillations presented in Fig. 
3. Such an intermittent behaviour can be presented by temporal frequency 
diagram [22,23] which can be given by the wavelet power spectra [24]. In Fig. 
6 one clearly see the laminar (periodic) interplay with turbulent (chaotic) 
regions. Note that for smaller rotation velocities we observe longer laminar 
intervals. 
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Fig. 6. Wavelet power spectra for torsional angular displacements (p; lo = 193, 200, 
201, 132 for (a), (b), (c), and (d); respectively. Note, the logarithmic colour scale 
on the right hand side. White regions in figures (a-d) indicate that power spectrum 
values exceed the scale. 



Table 1 
Parameters for a simple-supported Jeffcott rotor with a crack. 



Disk mass, M 


3 kg 


Disk polar moment of inertia, Jp 


0.01 kg m2 


Eccentricity of the disk, e 


2.2 X 10-5 m 


Torsional natural frequency, oot 


630 rad/s 


Lateral natural frequency, uot 


210 rad/s 


Shaft lateral stiffness, k 


k = LofM 


Shaft torsional stiffness, kt 


kt = uj^Jp 


Damping ratio in lateral direction, ^i (Ci = 2^iMuJi) 


0.008 


Damping ratio in torsional direction, ^t {Ct = '^^iJp^t) 


0.0006 


External torque amplitude, Tq 


800 Nm 


External torque amplitude, uj^ 


G [100, 250] rad/s 


Relative stiffness change caused by cracks — ^ 


0.4 


(Afci = |A%/fc, AA;2 = 1^) 





4 Discussion and Conclusions 



Using a simple model of the breathing crack we have showed that the response 
to periodic torque applied at the ends of rotating shaft is complex. In our re- 
sults we observe on one hand the multi-frequency synchronization, and on the 
other hand quasi-periodic motion and nonlinear beats in the regions of reso- 
nances. Especially the beats are characteristic as far as the nonlinear coupling 
of lateral and torsional modes induced by the crack is concerned. The beats are 
visible in both channels. They are also the signatures of interplay of external 
and parametric excitations, which arise from the crack existence. Interestingly 
the beats causing the intermittency of the laminar phases introduce the non- 
periodic element of vibrations. The chaotic nature of such dynamic response 
of the rotor system exceeds the present report and will be discussed in detail 
in a future publication. 
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